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Abstract 

The dynamics of electrons in the presence of a positive ion is considered for conditions of weak 
electron-electron couping but strong electron-ion coupling. The equilibrium electron density and 
electric field time correlation functions are evaluated for semi-classical conditions using a clas- 
sical statistical mechanics with a regularized electron-ion interaction for MD simulation. The 
theoretical analysis for the equilibrium state is obtained from the corresponding nonlinear Vlasov 
equation. Time correlation functions for the electrons are determined from the linearized Vlasov 
equation. The resulting electron dynamics is described in terms of a distribution of single electron- 
ion trajectories screened by an inhomogeneous electron gas dielectric function. The results are 
applied to calculation of the autocorrelation function for the electron electric field at the ion for 
< Z < 40, including conditions of strong electron-ion coupling. The electron stopping power 
and self-diffusion coefficient are determined from these results, and all properties calculated are 
compared with those obtained from semi-classical molecular dynamics simulation. The agreement 
with semi-classical MD simulation is found to be reasonable. The theoretical description provides 
an instructive interpretation for the strong electron-ion results. 

PACS numbers: 52.65.Yy, 52.25.Vy, 05.10.-a 



I. INTRODUCTION 



The total electric field at a particle in a plasma determines the dominant radiative and 
transport properties of that particle. The theory for the equilibrium distribution of fields at 
a neutral or charged point is well developed Q]. The theory for the dynamics of such fields 
is more complicated and some progress has been made recently in special cases. For 

example, the field dynamics at a neutral point has been described exactly in the Holtzmark 
limit j.2:] . The dynamical properties of fields due to positive ions near a positive impurity have 
been given an accurate approximate evaluation for a wide range of charge coupling, relative 
charge numbers, and relative masses 0, 0]. The corresponding study of negative charges 
(electrons) at a positive ion has been considered more recently only for the simplest case of a 
single ion of charge number Z in a semiclassical electron gas p,!^. In this case, the attractive 
interaction between the electron and ions emphasizes further the nonlinear dependence on 
Z. The static properties (electron charge density, electron microfield distribution) have 
been discussed in some detail for this case elsewhere Here, attention is focused on the 
dynamics via the equilibrium electron electric field autocorrelation function. This case of 
electron fields at a positive ion is qualitatively different from same sign ion fields, since in 
the former case electrons are attracted to the ion leading to strong electron-ion coupling for 
the enhanced close configurations. 

It is difficult a priori to predict even the qualitative features of the field autocorrelation 
functions due to this inherent strong electron-ion coupling, and there is no phenomenology 
for guidance. Consequently, our initial analysis has been based on MD simulation of the 
correlation functions followed by an attempt to model and interpret the observed results. 
However, MD simulation for the electrons is limited to classical mechanics, while the singular 
attractive electron-ion interaction inherently requires a quantum mechanical description. 
This difficulty is circumvented by modifying the electron-ion Coulomb potential at short 
distances to represent quantum diffraction effects. The conditions for validity and limitations 
of this classical model have been discussed extensively elsewhere 9| . The details of the MD 

Q 

method also have been discussed elsewhere j7|] and will not be repeated here. There is a 
growing recent literature on the related MD studies of two-component classical models of 



a hydrogen plasma 
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ll| at stronger electron-electron coupling values, but restricted to 



Z = 1. Thus the main new feature studied here is the dependence of structure and dynamics 



2 



on charge number Z. 

The relevant dimensionless parameters are the charge number of the ion, Z, the electron- 
electron coupling constant F = e^/rQkBT, and the de Broglie wavelength relative to the 
interelectron distance, 5 = {2Txh^ / rrif.kBTr'lY/'^ . The interelectron distance is defined in 
terms of the electron density rif. by 47rnerQ/3 = 1. The electron- ion coupling is measured by 
the maximum value of the magnitude of the regularized electron-ion potential at the origin, 
a = Zr/6. Most of the results described below are for F = 0.1, 6 = 0.4, and a = 0.25Z, 
with < Z < 40. The corresponding density and temperature are = 2.5 x 10^^ cm~'^ and 
T = 7.9 X 10^ K. Additional results are presented for the (unrealistic) extreme electron-ion 
coupling conditions of F = 0.5, 6 = 0.2, and a = 2.5Z, with < Z < 10 (n = 3.2 x 10^^ 
cm~^, T = 7.9 X 10^ K ), and for the experimentally relevant conditions of F = 0.029, 
S = 0.059, and a = 0.5Z, with 0<Z<8(n = lx 10^^ cm'^, T = 2 x lO^K ). In 
all cases the electron-electron coupling is weak. Since F is small the kinetic equation for 
the electron reduced distribution function becomes the nonlinear Vlasov equation in the 
presence of the external ion potential. For the equilibrium state this equation gives the 
nonlinear Boltzmann-Poisson equation. It is mathematically equivalent to the HNC integral 
equations for an impurity in a one component plasma applicable as well for larger F, and 
can be solved numerically in a similar way. However, as noted elsewhere ||7|, these numerical 
methods fail at very strong coupling (cr > 8). For practical purposes, the stationary solution 
is modelled as a nonlinear Debye-Huckel distribution with parameters fit to the HNC solution 
when it exists or by comparison with MD simulations otherwise (a variational method also 
has been described but is not used here). The results for equilibrium properties are 
shown to be quite accurate. 

The equilibrium time correlation functions are determined from the linearized Vlasov 
equation (linear perturbations of the equilibrium state, but still nonlinear in the electron- 
ion interaction). The results are a composition of correlated initial conditions, fields with 
single electron trajectories about the ion, and dynamical collective screening by the inhomo- 
geneous electron distribution about the ion. These results are simplified further by a mean 
field approximation resulting in a single electron problem in the effective potential of the 
nonlinear Debye distribution. This analysis is applied to the electric field autocorrelation 
function showing reasonable agreement with the results from MD simulation. The primary 
observations from the simulations of the field autocorrelation function for increasing charge 



3 



number are: 1) an increase in the initial value, 2) a decrease in the correlation time, and 
3) an increasing anticorrelation at longer times. The simple mean field approximation re- 
produces all of these results. Furthermore, its simplicity allows an interpretation from the 
single particle motion. 

A closely related property is the stopping power for an ion by the electron gas. In the 
low velocity limit, and for large ion mass, the stopping power is proportional to the time 



integral of the field autocorrelation function 



li 



diffusion and friction coefficients in this same limit 



'his integral also determines the self- 



la] . Linear response predicts a dominant 



dependence for these properties. Significant deviations from this dependence are 
observed at strong coupling and have been the focus of attention in recent years 2. The 
results here show these deviations come from a competition between the increase of the 
integral due to 1) above and the decrease due to 2) and 3). MD simulation shows that the 
latter two dynamical effects dominate the former static effect. The simple mean field model 
provides the missing interpretation for 2) and 3). 

This same integral of the field autocorrelation function determines the half width for 
spectal lines from ion radiators due to perturbations by electrons in the fast fiuctuation 



limit 



21[. An extension of the analysis provided here to this atomic physics problem is 



under way and promises to provide an additional experimental probe for the dynamics of 



charges near an ion 
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The theoretical description based on the Vlasov equation is provided in the next section. 
The results are applied to the field autocorrelation function in section 3. The stopping power 
and self-diffusion coefficient are evaluated in section 4. Finally, a summary and discussion 
is provided in the last section. 

II. KINETIC THEORY 

The classical system considered consists of A'^e electrons with charge — e, an infinitely 
massive positive ion with charge Ze placed at the origin, and a rigid uniform positive back- 
ground for overall charge neutrality contained in a large volume V. The Hamiltonian has 
the form 

^ = E ( ^ + (^") + (^-)] + 2 5Z ^-d^" - (1) 
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where Ke(|ra — r^|)is the Coulomb interaction for electrons a and 7, Veiira) is the reg- 
ularized electron-ion interaction for electron a with the ion, and Veb is the Coulomb 
interaction for electron a with the uniform neutralizing background 



e 



2 



(1 - e-''"/''0'5) 



Vee{\ra ' r^|) = , _ K^(r«) = -Ze' ^ '-. (2) 



V;6(ra) = -(^-^^n, j drKedra - r|). (3) 

For values of r/ro >> 6 the potential Vei{r) becomes Coulomb, while for r << 6 the Coulomb 
singularity is removed and (3V^^{r) —f3Ze^/rQ5 = — = —ZT/5 = —a. This is the simplest 
phenomenological form representing the short range effects of the uncertainty principle Q| . 
In principle there should be a similar regularization of the electron-electron interaction, but 
since that interaction is repulsive configurations with a pair of electrons within a thermal 
de Broglie wavelength are rare. For simplicity, therefore, the electron-electron interaction 
is taken to be Coulomb. In all of the following the electron-electron coupling (weak) is 
measured by F while that for the electron-ion coupling (possibly strong) is measured by a. 

The distribution function for the electrons is denoted by /(r, v; t). It is normalized such 
that integration over all velocities and the volume of the system equals N^,. It obeys the 
exact first BBGKY (Born, Bogoliubov, Green, Kirkwood, Yvon) hierarchy equation [l7| 

(9i + V ■ Vr - m-^ (Vr (K.(r) + K6(r))) ■ Vv) /(r, v; t) 
= m-^ j rfr2rfv2(VrV;e(r-r2))- Vv/('Hr,v;r2,V2;t) (4) 

Here /^^-"(r, v; r2, V2; t) is the joint distribution function for two electrons. At the 
weak coupling (F << 1) the electron distributions are approximately independent and 
/'^^)(r, v; r2, V2; t) /(r, v; t)/(r2, V2; t). Then ()H) becomes the nonlinear Vlasov equation 

(a* + V ■ Vr - m-^ (Vr (K.(r))) • Vv) /(r, v; t) 
= m-' (Vv/(r, v; t)) ■ Vr j dr^VUr - r2)5n(r2, t) (5) 

where (5n(r, t) is the deviation of the electron density from the rigid uniform positive back- 
ground 

6n{r, t) = n{r, t) - (^1 - n^, n(r, t) = j dvfi^, v; t). (6) 

The left side of (0) describes the single electron motion in the presence of the ion at the 
origin. Since the charge number can be large, the electron ion coupling can be large. The 



right side of this equation describes the correlation effects due to interactions among the 
electrons, which are weak for hot, dense matter. However, these weak correlations depend 
nonlinearly on the distribution function so that the content of this equation and its solutions 
can be quite rich and complex. 



A. Equilibrium solution 

The equilibrium state is a stationary solution to (0). However, it is known from classical 
equilibrium statistical mechanics that it must be of the form 

/,(r, v) = n,{r)cP{v), <Piv) = (^^^ e-^"''V2. (7) 

Substitution of this into (0) gives the nonlinear integral equation for ne{r) 

l^'lM ^ .pV„(r) - ^e^ [ dr'^J^ (8) 

He J \r — r I 

where rig = Ne/V is the average volume. It is often useful to write the solution in terms of 
an effective electron-ion potential 

He 

According to (jH)) Uei{r) obeys the nonlinear integral equation 



\r — r\ 

The second term provides the nonlinear strong coupling effects of the electron-ion interac- 
tions. It is worth noting that although the electron-electron coupling is weak, the strong 
ion-electron effects are mediated by the electron-electron interaction of this second term. 
The predictions for a free electron gas interacting with an ion are quite different from those 
given below. 

The numerical solutions to (jHl) or (jlO|) and its relationship to the hypernetted chain 
integral equation at stronger electron coupling has been discussed in reference ^\. At the 
weak electron-electron coupling considered here these equations are the same as the HNC 
equations, and their solution in the following will be referred to as the HNC result. An 
important observation in reference |3] is that the numerical solution is well represented by 
the Debye form 



-Ze_^ 1 

(5Af 



Ue^{r') = 7 e-'-/^«^ - e-'/''^') . (11) 
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FIG. 1: Electron density around an ion of charge number Z = 8 



where the effective charge number Z and screening length A are fitting parameters. In the 
weak couphng domain these become the actual charge number and the Debye screening 
length and this form is exact. More generally, it is only an approximation and the best 
choices for Z, X are different from Z, A at strong coupling. This approximation for ne{r) 
will be referred to as the nonlinear Debye approximation. Figure 1 illustrates the fitting of 
the nonlinear Debye form to HNC for the case Z = 8, F = 0.1, and S = 0.4. Also shown are 
the corresponding results from MD simulation. The agreement is quite good. Finally, the 
linear Debye form is shown to indicate that nonlinear effects are clearly significant. Similar 
results are obtained for the other values of Z discussed below. 
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B. Electric field time correlation function 



To explore the dynamics of the electrons attention is limited to the electric field autocor- 
relation function defined by 



C{t) = < E (t) -E > . 



(12) 



The brackets denote an equilibrium Gibbs ensemble average and E is the total field at the 
ion due to all electrons 

Ne Ne 

(13) 



E 



Ne Ne 

5Z = 5Z ^ '^^^ 



a=l 



e(r„) = e^(l-(l + ^)e--/-^) 



Consider first the initial value C(0) for which two exact representations can be given 

ne(r)e(r)+ /" c/r'rigfr, r')e (r') 

drnp(r)e (r) -e^ 



C{0) = ^ y rfre(r) 



(15) 



The first equality expresses the covariance in terms of both the electron density ne(r) and 
the density for two electrons near the ion ne(r, r'). The second representation requires only 
the electron density and the mean force field derived from it 



rieir) 

1 



-Vln 



ne(r)e (r) + J dr'ne{r,r')e (r') 

ne{r) 



(16) 



f3Ze ne 

The second equality is derived in Appendix A. It is similar to the first form with the apparent 
neglect of the two-electron ion correlations. However, these latter contributions are incorpo- 
rated exactly in the mean force field e^/ (r) . The effects of electron-electron interactions are 
included in ne(r) through the mean field screening of the ion-electron interaction. However, 
ne(r, r') describes the additional electron correlations for two electrons near the ion beyond 
these mean field effects. Figure 2 shows the initial field fiuctuations C(0) as a function of Z, 
at r = 0.1 and a = 0.25Z, calculated using ne(r) from HNC and calculated directly from 
MD. The solid curve is the HNC data indicating a strong ~ Z^ dependence inherited from 
the electron charge density. Also shown are the results from the first equality of p5|) neglect- 
ing the correlations of two electrons in the presence of the ion, ne(r, r') —>■ ne(r)ne(r')/ne. 
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Clearly, these are significant for the conditions considered, showing the coupling between 
strong ion-electron interactions and weak electron-electron interactions. The results are 
similar for the other coupling conditions considered (cr = 0.5Z and a = 2.5Z). 

It is shown in Appendix B that the time dependence of C{t) in the weak coupling limit 
can be written as 

C(t) = y"drdve(r) ■^(r,v,t) (17) 
where ^/'(r, v, t) obeys the linearized Vlasov equation 

(dt + C) ^(r, v; t) = -/5/e(r, v)v ■ Vr j cir2\4e(r - ra) j rfv2^(r2, V2, t). (18) 

The associated initial condition is 

^(r,v,t = 0) = e™/(r)/e(r,v), (19) 

where em/(r) is the mean field of (fTBj) above. The linear operator C is the generator for 
single electron dynamics in the effective potential Uei{r) 

£ = V • Vr - m;^ (Vr {Ue^{r))) ■ Vv (20) 
9 



The equilibrium distribution is stationary under this operator, £/e(r, v) = 0. In addition to 
this effective single particle dynamics, all dynamical many-electron effects are contained in 
the term on the right side of (|18|). The details of the formal solution to this kinetic equation 
for the correlation function are given in Appendix B with the result 

^^t) = j^^' J drd^rf,{r,^r)e,{r■t')e^Jirit -t')) (21) 

The field e(r(t — t')) is given by ()14|) with the initial position shifted to r{t — t') according to 
the effective single particle dynamics generated by C, using the initial conditions r, v. The 
other field es(r; t') is a dynamically screened field 

e,(r;t) = C^^)'^ j dke-'^-^e,{k,t), e,(k,t) = j dk'e{k')e-^ (k',k,t) (22) 

where e(k') is the Fourier transform of (fT^ . and the dielectric function is defined by 

e(k, k', t) = {2nf6(k - k') + 7r(k, k', t)Vee{k') (23) 



n{k,k',t) = /" rfrrfv/e(r, v)e*''"e-*^'-'"W (24) 

ut J 

For Z = this is the familiar classical random phase approximation for a uniform electron 
gas, diagonal in k, k'. More generally, the Z dependence leads to a nonuniform electron 
density near the ion and the polarization function 7T(k,k',t) depends on the details of this 
distribution. It vanishes at t = indicating no initial screening, 7r(k, k', t = 0) = 0. At later 
times the polarization function is non-zero giving a space and time dependent additional 
screening. Further simplification of this result for practical evaluation is discussed below. 



III. MOLECULAR DYNAMICS SIMULATION 



The application of standard MD simulation methods using the semi-classical electron-ion 
potential is somewhat more complex than for the usual classical fluids with short range 
repulsive interactions. Although finite at short distances the attractive electron-ion poten- 
tial allows bound and metastable states for electrons orbiting round the ion over extended 
periods. For most properties, e.g., structural properties, this is not a severe problem except 
at low temperatures. However, the interest here is in electric field dynamics which is very 
sensitive to close electron-ion configurations. The protocol for control of anomalous states 
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FIG. 3: Field autocorrelation function at strong coupling 

with quasi-bound dynamical states has been described in reference j3] and will not be re- 
peated here. Consider first the coupling conditions T = 0.1 and a = 0.25Z. Figure 3 shows 
the field correlation function C{t) for Z = 4,8, 20, 30, 40. There are two qualitative features 
to note. The first is a characteristic time for relaxation that decreases with increasing Z, 
and the second is the development of an anti-correlation that increases with increasing Z. 
Figures 4 and 5 show the corresponding results for the weaker ( F = 0.029 and cr = 0.5Z, 
for Z = 1,3,5,8) and stronger ( F = 0.5 and cr = 2.5Z, for Z = 1,2,3,4,6,10) coupling 
cases. In the former case the decreasing relaxation time is evident but the anti-correlation 
is significant only for the Z = 8 curve. The strongest coupling case of Fig. 5 shows large 
anti-correlation for all Z. It should be noted that this last case is unrealistic since the 
equilibrium population of ions at such strong coupling is extremely small. 

The area under the curves is related to transport properties, as indicated in the next 
section, and results from a competition between these two features and the increasing initial 
correlations. The effects of shortening decay time and anti-correlation dominate to decrease 
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the area as indicated in Fig. 6. These effects are greater for stronger couphng, with sig- 
nificant anti-correlation occurring for Z > 4. Figure 7 shows the same results plotted as a 
function of a. The three cases appear similar, with only a shift in their a = value due to 
different values for 7 and 5. This suggests that the time integral of the correlation function, 
normalized to its value at cr = may be a "universal" function of sigma. 

IV. EFFECTIVE SINGLE PARTICLE MODEL 

To calculate the correlation function and identify the mechanisms for the decay time 
and correlation, it might be supposed that the single nearest electron dominates since its 
field is greatest. Figure 8 shows this nearest neighbor contribution from MD for the same 
conditions of Fig. 3. Clearly, the nearest neighbor field autocorrelation function has the same 
qualitative behavior as the total correlation function with respect to the decreasing decay 
time and increasing anti-correlation. The quantitative values are wrong (both amplitude and 
time scale), however, suggesting that contributions from other particles are important as well. 

12 
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FIG. 5: Field autocorreltion function at very strong coupling 



To explore this in more detail, the Vlasov equation given by pij) can be used. It requires 
evaluation of the dynamically screened field es{r;t), although at short times es{r;t) ~ 
e(r; 0) = e(r). The correlation function is then effectively that for a single particle moving 
in the self-consistent potential (fTIUl . averaged over the initial equilibrium distribution of 
electrons about the ion. Thus it is similar to the nearest neighbor approximation but extends 
it to include all electrons, including the correct initial correlations and the correlations of 
the mean field for the dynamics. If the additional dynamical screening is neglected for all 
relevant times, i.e. 

e,(r;t) ^5(t)e(r) (25) 

then (PT|) becomes 

C{t) ^ J rfrrfv/,(r,v)e(r)e^^(r(t)). (26) 

Comparison with (fT^ shows that this approximation is exact for C(0). For practical pur- 
posed, this approximation for C{t) has been evaluated using the analytic nonlinear Debye 
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FIG. 6: Time integral of C{t) as a function of Z for moderate, strong, and very strong coupling 
cases 

form (fTT|) fitted to the HNC results in both the dynamics and the electron density. The 
results for this effective single particle model are shown in Figs. 9 and 10 for the case F = 0.1 
and a = 0.25Z. 

The agreement between this simple theory is quite good and provides a means to interpret 
the results of the MD simulation. The initial position and velocity of an electron are sampled 
from the equilibrium distribution /e(r, v) which favors electrons close to the ion and hence 
large fields. Since the force on the electron is also large, its initial acceleration will be large. 
This is the source of the short decay time. Consider the case of an energetic electron near 
the ion. If its velocity is directed away from the ion it will move to larger distances and 
the field will decrease with a positive component along the initial field. In contrast, if its 
initial velocity is toward the ion it will move past the ion with a change in the direction of 
its field relative to the initial value. This is a source for anti-correlation. For less energetic 
electrons, the trajectories are bound and there is continual correlation and anti-correlation 
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FIG. 7: Same as Fig. 6 function of a 



as the correlation function decays in magnitude due to phase averaging. Botli the increase 
in initial correlation and the field reversal effects should increase as the charge on the ion 
increases, and this is what is observed. At weaker electron-ion coupling both effects are 
diminished and blurred as the relevant configurations are more distant, the fields are weaker, 
and the accelerations smaller. This is already evident in the nearest neighbor results of Fig. 
8. 

The simple model of (I26p appears better at larger values of Z where the single particle 
motion is expected to dominate. For smaller values of Z the agreement at short times is 
still good (the discrepancy at t = is a limitation of HNC, not the dynamics), but more 
significant differences occur after the first initial decrease. Presumably, this is due to the 
dynamical screening effects in es(r;t) that have been neglected. It should be noted that 
the results here are somewhat sensitive to the choice of parameters Z, A used in fitting the 
non- linear Debye Huckel form for /e(r, v). Fits emphasizing short or intermediate distances 
change slightly the point at which anti-correlation sets in and its amplitude. The primary 
criterion used here was a globally good visual fit and a good resulting value for the initial 
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FIG. 8: Nearest neighbor field autocorrelation function for strong coupling 



condition C(0). 

V. STOPPING POWER, FRICTION, AND SELF-DIFFUSION 

Emphasis here has been placed on the electric field autocorrelation function as a sensitive 
measure of electron properties near the ion. This function is also of interest because of 
its connection to transport and radiative properties of the ion. Specifically, for the case 
of an infinitely massive ion considered here there are exact relationships between transport 
coefficients characterizing three physically different phenomena: 1) the low velocity stopping 
power S for a particle injected in the electron gas, 2) the friction coefficient C, for the 
resistance to a particle being pulled through the gas, and 3) the self-diffusion coefficient D 
of a particle at equilibrium with the gas [15[ 

= iPny' = ^ |,=o= ^ r dtC{t) (27) 

^ Jo 
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FIG. 9: Comparison of C{t) from MD with results from mean field kinetic theory for Z = 4,8 



Finally, the time integral of C{t) also provides the fast fluctuation limit (impact) for the 
spectral line width of ions broadened by electrons Clearly, a better understanding of 
the mechanisms controlling the electric field autocorrelation function is of interest in several 
different contexts. 

This Green-Kubo representation ()27p allows a determination of these transport properties 
from an equilibrium MD simulation, as described above. In contrast, previous simulations 
of stopping power have studied the nonequilibrium state of the injected particle, measuring 
directly the energy degradation |I9] . At asymptotically weak coupling, these properties have 
a dominant dependence, as C{t) becomes independent of Z. A puzzling result of the 
previous simulations and some experiments j^, was the observation of a weaker Z 
dependence at strong coupling. This behavior is somewhat puzzling in light of the strong 
growth of the initial value C(0) ~ Z^ at large Z (see Fig. 2). Thus it would appear that 
the dominant dependence would be an even stronger Z^. However, as Fig. 6 shows clearly 
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FIG. 10: Same as Figure 9 for Z = 20, 30 

the competing effects of decreasing correlation time and a developing time interval of anti- 
correlation dominate at strong coupling to decrease the time integral. Figure 11 shows the 
dimensionless stopping power as a function of Z for the case F = 0.1 and a = 0.25Z. Also 
shown is the Born approximation 0.83Z^, where the coefficient has been determined by the 
data for small Z. The MD data has been fit to a crossover function 

0.83^2 



Siv) 



\v=0' 



(28) 



1 + 0.008^1/2- 

This form has been chosen since it implies the stopping power goes as Z^^"^ at extreme 



coupling, which is consistent with the earlier results [19|, |20|. However, other fits to the data 
here are possible as well. The predictions of the simple effective single particle theory are 
shown on Fig. 11 also. 

Similar results are obtained for the other coupling cases. Figure 12 shows the stopping 
power for the weaker coupling case of F = 0.029 and cr = 0.5Z at smaller values of Z. 
The Born approximation l.SSZ"^ is determined by the Z = 1 value. As expected the Born 
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FIG. 11: Stopping power relative to the speed at zero speed, Eq. (|27j) . at strong coupling. Also 
shown is the Born approximation and the results of the kinetic theory 



approximation is quite good at weak coupling, although some deviation is seen at Z = 8. 
This is consistent with the above estimate that strong coupling effects in Fig. 6 occur for 
(T > 4. The strongest coupling case of F = 0.5 and a = 2.5Z is shown in Fig. 13, where the 
Born approximation is again determined from the Z = 1 data and found to be 1.52Z^. The 
deviations from the Born approximation are much greater now, as expected. 



VI. SUMMARY AND DISCUSSION 



The objective here has been to explore the dynamics of electrons near a positive ion as a 
function of the charge number on the ion, or more precisely, as a function of the electron-ion 
coupling a = ZT /5. A primary tool for this investigation has been molecular dynamics sim- 
ulation, requiring a semi-classical regularization of the Coulomb potential at short distances. 
Under the hot, dense conditions needed to support large Z ions the electron coupling can 
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FIG. 12: Same as Fig. 11 at moderate coupling 



be quite weak. An accurate theoretical description is then given by the nonhnear Vlasov ki- 
netic equation for the electrons. The MD simulation reveals an interesting structure for the 
electric field autocorrelation function. At the weakest coupling considered there is a rapid 
initial decay followed by an asymptotic domain of weak anticorrelation. The time scale for 
the initial decay decreases with increasing Z but does not depend strongly on Z as long as 
the coupling stays weak. This is illustrated in Fig. 4. Figures 5 and 6 show a quite different 
behavior at strong coupling where the initial decay time shortens and the anticorrelation 
becomes prominent on the same time scale. This qualitative behavior is characteristic of 
the single electron dynamics of the nearest neighbor. A quantitative description is provided 
by the Vlasov kinetic theory with the exact initial condition for the distribution of electrons 
about the ion, including strong electron-ion correlations mediated by weak electron-electron 
interactions. The mean field dynamics is that of effective single electron trajectories, calcu- 
lated for the same effective potential as that for the initial correlations, and averaged over 
an ensemble of these initial states. 
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FIG. 13: Same as Fig. 11 at very strong coupling 



The simple theoretical model provides a means to interpret the MD simulation data for 
the beginnings of a phenomenological understanding of electron dynamics near a positive 
ion. The initial correlation for the electric field C(0) increases with Z as the equilibrium 
distribution of electrons is enhanced near the ion with configurations corresponding to larger 
fields. The latter is well described by the nonlinear Debye-Huckel form given by (jH)) with 
(jlll) . The initial decay of Of) is essentially the decorrelation time for a "most probable" 
electron near the ion. This most probable distance can be estimated from the maximum 
of the Debye distribution P(r) = 47rr^ne(r)/ne = 47rr^e~^'^"^^^ to give r/rg ~ b for strong 
coupling. The correlation time is then approximately the time to accelerate this electron to 
the position of the ion, t^. ~ (ro^/i/^m) ^J\/ZT. The first factor r^b/ \/(3m is the time for 
an electron with the thermal velocity to cross a sphere of the size of the thermal de Broglie 
diameter. The second factor gives the dominant Z dependence tc ~ ^Jl/Z . As the electron 
continues past the ion its acceleration changes sign as it is attracted to the ion with a field 
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opposite that of the original field. This is the source of the dominant anticorrelation. 

The self-diffusion coefficient, stopping power, friction coefficient, and width of spectral 
lines are proportional to the integral of C(t). At Z = this is the field autocorrelation 
function at a neutral point and is independent of Z, depending only on the electron-electron 
coupling. As Z increases the electron distribution becomes nonuniform about the ion and 
C(0) increases. At the same time tc decreases. For small Z these two effects become negligi- 
ble, as seen in Fig. 6 for the case of cr = 0.5Z and Z < 5. As Z increases, or more precisely 
as a increases, the domain of anticorrelation appears and begins to dominate the decrease 
in the integral of Of). These basic mechanisms are captured by the meanfield description 
based on the Vlasov equation for the one particle electron distribution. The relevant corre- 
lations contained in this description are those of the equilibrium electron distribution about 
the ion, also described well by the stationary solution to the Vlasov equation. 

This analysis provides a new picture for the puzzling decrease of stopping power with 
increasing Z, relative to the Born approximation. The stopping power is proportional to Z^ 
times the time integral of C(t) which is essentially the total cross section for all the electrons 
and the ion. The decrease towards a Z"^!"^ dependence at strong coupling observed earlier is 
seen to be due to the effects just described. However, the precise dependence on Z may be 
more complicated as the coupling increases. 

Similarly, these same results provide clear evidence for the effects of electron-electron 
and electron- ion correlations on the shape of spectral lines 22 , [2^ . Similar experimental 



puzzles regarding the Z dependence of the half width |2J| can be clarified through combined 
theoretical and MD simulation as described here. 
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APPENDIX A: EVALUATION OF THE FIELD COVARIANCE 



The dimensionless electric field covariance is defined by 



C(0) 



< E ■ E > 



with 



E = ^ e (r« - ro) . 



(Al) 



(A2) 



a=l 



Here Yq denotes the position of the ion. The equilibrium average can be calculated directly 
in terms of the one and two electron charge densities 



C(0) = ^ j rfre(r)- ne(r)e(r) + j rfr'ne(r, r')e (r) 



with the definitions 



where ?7 is the total kinetic energy and T = l/ksP is the temperature. 
An equivalent alternative form is obtained by writing the covariance as 



(A3) 



(A4) 



C(0) 



<VroUie■E>-- 



J dvodv2..drN (V.^e-^^) -E 



J dvQ.AvNe' 



■pu 







(Vn, ■ E) 



drnp(r)V -e (r) 



J drueir) l^-^^Vlnrae(r)^ •e(r) = ^ J rfrne(r)e„/ (r) -e (r) (A5) 



In the last equality the mean force field has been introduced by 

1 



pZe 



V In ne(r, t). 



Finally, comparison of ()A3|1 and ()A5|1 gives the alternative expression for this field 

1 



ne{r) 



ne{r)e{r) + J (ir'ne(r, r')e (r') 



(A6) 



(A7) 



APPENDIX B: KINETIC EQUATION FOR CORRELATION FUNCTIONS 

In this section the evaluation of the field autocorrelation function by kinetic theory, and 
the basis for the approximation ()18|). are briefly described. First, the correlation function is 
formally rewritten as 
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4 4 /• 

C{t) = ^<E(t)-E>=^ / dridvi..cir^dvjvE- E(-t)pe 
-| J dridvie{ri) N- J dr2dv2..dvj^dvNE (—t) pe 



e 

A 



4 /* 

= ^ / (iri(ivie(ri)-^(ri,vi;t) (Bl) 

where pe is the equihbrium Gibbs ensemble and e (r^) is the single particle field of p8|) . The 
integrations over degrees of freedom L.N in the second equality define a reduced function 
■^(ri, Vi; t) which is the first member of a set of such functions 

V^(")(ri,vi; ..r„ v,;t) = A^" ^ c/r,+ic/v,+i..c/rjvc/v^E (-t) pe- (B2) 

It is straightforward to verify that these functions satisfy the BBGKY hierarchy, whose first 
equation is formally the same as (@]) 

(5< + V ■ Vr - (Vr (K.(r) + Kblr))) ■ Vv) ^(ri, Vi; t) 
= m;' j rfr2rfv2(VrV;e(r-r2))- VvV^('nr,v;r2,V2;t). (B3) 

However, in contrast to the distribution functions in (@]) the functional relationship of ip'^'^^ 
to ijj is linear. To see this, consider first the initial conditions which are found to be 

V'(ri,vi;t = 0) = /e(ri,vi)e(ri) + j (ir2(iv2/i^^(ri, Vi; ra, V2)e (r2) 

= /e(ri,vi)e„y(ri) (B4) 

^(2)(ri,vi;r2,V2;t = 0) = /f (n, Vi; r2, V2) (e (n) + e (r2)) 

+n J (ir3(iv3/f)(ri, vi;r2, V2;r3,V3)e(r3) 

= /i^nri,vi;r2, V2) (e„ (ri) + e„ (r2)) 

+n j (ir3dv3/i^^)(ri,vi;r2, V2;r3, V3)e(r3) (B5) 

Here, /i**"* are the equilibrium s— particle reduced distribution functions associated with the 

(3) 

Gibbs ensemble and he is the equilibrium correlation function for three electrons in the 
presence of the ion 

/if)(ri,vi;r2,V2;r3, V3;t) = /f ^(ri, Vi; r2, V2; r3, V3) 

-/i^^ (ri, vi; r2, ^2)^, r/i^^ri, Vi; rg, V3) 

/e(ri,vi) 

-/f (ri,v,;r2,V2)— (r2,V2;r3,V3) (B6) 

Je[r2, V2j 
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The linear functional relationship between ip^"^^ to at t = is now evident 

/(2)/ n^ /(2)/ ,x /^(ri, vi;0) ^(r2,V2;0) 
^(ri,vi;r2,V2;0) = ^(ri, vi; r2, V2; t) - + 



/e(ri,Vi) /e(r2,V2) 

+3 electron correlations. (B7) 

Recognizing this linear relationship, the basic approximation for weak coupling among 
the electrons is to neglect all of their correlations at all times, i.e. extend flB7|) to 

V'^^Hri, vi;r2, V2;t) /e(r2, V2; t)?A(ri, Vi; t) + /e(ri, Vi; t)V^(r2, V2; t). (B8) 

Use of this in the first hierarchy equation ()B3|) gives directly the kinetic equation ()18|) 
discussed in the text 

{dt + C) ^Pir, v; t) = -/3/e(r, v)v ■ Vr J dr2Veeir - r^) j rfv2^(r2, V2, t). (B9) 



£ = V • Vr - m-' (Vr {U,,{v))) ■ V^. (BIO) 



APPENDIX C: SOLUTION TO KINETIC EQUATION 



The operator C in ()B3jl is the generator for single electron dynamics in the effective 
potential due to the ion Uei{r). The solution to the equation can be obtained in terms of 
this single electron dynamics by direct integration 

^(r, v; t) = e-^V(r, 0) - ^ rfre-^(*-^)/e(r, v)/?v ■ Vr j dv^Vee^v - r2)I(r2, r) (CI) 

I(r,t) = j rfvV^(r,v,t) (C2) 
The initial condition is given by ()B4|) . An equation for I(r,t) follows from substitution of 

dnu into (ini 

I(r;t) = j rfve-^Ve(r,v)e^^(r)- j' dr j rfve-^(*-^)/e(r, v)/5vVr ^ c?r2V;e(r - r2)I(r2, r). 

(C3) 

This is an integral equation for I(r; t) which can be written 

dr j dr2e{r,t;r2,T)I{r2;T) = y c/ve"^*/e(r, v)e™j(r). (C4) 
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The dielectric function e (r, t — r; is defined by 

e (r, t - r; r') = 5 (t - r) 5 (r - r') + j dr"n (r, t - r; r") Veeir"-r'), (C5) 
where the polarization function is 

TT (r, t; r') = J rfve-^*/e(r, v)/3v ■ Vr5 (r - r') . (C6) 
With these results the correlation function from ()B1|) becomes 
= ^ / drdve{r)-ij{r,v;t) = ^ / c/re (r) ■ I(r; t) 

= j ^r^^*^^ ^) ■ e-'^^^-^Velr, v)e„,Xr). (C7) 

In the last equality the screened field (r, t) has been introduced 

e,(r,t) = j dr'e(r')e-i(r',t;r) (C8) 
Finally, using the stationarity of /e(r, v) under the dynamics generated by C gives 
C{t) = "1 j'^dr j rfrrfv/e(r,v)eJr,r)-e-^(*-^)e„Xr) 

This is the form used in the text, Eq. (j2H) . 
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